Arithmetic means of speculative coverage are in Figure 1.

Figure 1: Arithmetic means of crop speculative coverage.

coverage.lm1 <- lm(speculative.pct.coverage   ~ as.factor(block) + species , data = coverage, subset = (treatment != "N"))
resid_panel(coverage.lm1)

joint_tests(coverage.lm1)

##  model term df1 df2 F.ratio p.value
##  block        3  30   8.800  0.0002
##  species     10  30  12.567  <.0001

plot(emmeans(coverage.lm1 , "species"), comparisons = TRUE) + 
  xlab("Speculative coverage (%)" ) +
  ylab("Species")  + 
  theme(text = element_text(size=20))

Figure 2: Estimate speculative percent coverage of different cover crop species